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Abstract 

The Toda lattice hierarchy with self- consistent sources and their Lax representation 
are derived. We construct a forward Darboux transformation (FDT) with arbitrary 
functions of time and a generalized forward Darboux transformation (GFDT) for Toda 
lattice with self-consistent sources (TLSCS), which can serve as a non-auto-Backlund 
transformation between TLSCS with different degrees of sources. With the help of 
such DT, we can construct many type of solutions to TLSCS, such as rational solution, 
solitons, positons, negetons, and soliton-positons, soliton-negatons, positon-negatons 
etc., and study properties and interactions of these solutions. 
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1 Introduction 



Recently, soliton equations with self-consistent sources (SESCSs) in several 1+1 and 2+1 
dimensional continuous cases, which have very important applications in many fields, such 
as hydrodynamics, solid state physics, plasma physics, etc. have been widely studied [1-5]. 
For example the KdV equation with self-consistent sources describes the interaction of long 
and short capillary-gravity waves [4]. Various methods have been used to construct their 
solutions, such as Inverse scattering methods [2,3,6,7], Darboux transformation methods 
[8-11], Hirota bilinear methods [12-14] etc. Comparing with Darboux transformations (DT) 
for soliton equations, generalized binary Darboux transformations (GBDT) with arbitrary 
functions of time developed in [8-11] provides a non-auto-Backlund transformations between 
two soliton equations with different degrees of sources, and enable us to find various solutions, 
such as soliton, positon, negaton, soliton-positon, soliton-negaton, positon-negaton and etc. 

However, in contrast with the continuous case, the integrable discrete systems with self- 
consistent sources and their physical applications have not been studied yet. We will present 
a way to construct integrable discrete systems with self-consistent sources basing on the con- 
strained flows of the integrable discrete system [15] by regarding the latter as the stationary 
equation of the former. Then we derive a forward Darboux transformation (FDT) with an 
arbitrary function of time for the integrable discrete system with self-consistent sources by 
improving the method in [8-11]. We use a discrete system, the Toda lattice equation with 
self-consistent sources (TLSCS) to illustrate our method. The FDT with an arbitrary func- 
tion of time is generated by using a linear combination of two independent eigenfunctions 
of the Lax pair with the combination coefficient explicitly depends on time. It serves as 
a non-auto-Backlund transformation between two TLSCS with different degrees of sources. 
We give a formula for multi-time repeating of such DT. We also construct the generalized 
FDT (GFDT) with arbitrary functions of time. The formula for GFDT is quite similar to 
the generalized Wronski determinant formulae [16] and the generalized Casorati determinant 
formulae [17,18]. 

It is well known that Toda lattice equation possesses rich families of solutions including 
rational solution, soliton, positon, negatons and soliton-positon, soliton-negaton, positon- 
negaton. Soliton is a fast decaying pulse like solution without singularity [21]. It describes a 
wave propagating in a constant speed. Positon is an oscillating and slowly decaying solution 
with singularities. It leads to a trivial scattering matrix (called " supertransparent" ) when 
inserted as potentials in the finite-difference Schrodinger equation of the corresponding Lax 
pair [17]. The interaction between positon and other types of solutions is very interesting. 
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There is no phase shift for others during the course of collision. Negatons is another type 
of solutions with singularities oscillating but fast decaying and with non-trivial scattering 
matrix. Unlike the positon case, there is phase shift for other types of solutions during the 
interaction between negaton and others. 

In the second part of our paper, we show that TLSCS also have such types of solutions. 
These solutions can be obtained easily through GFDT. Some of these solutions as well as their 
analytic properties are presented in our paper. Differences between solutions for TLSCS and 
solutions for Toda lattice equation are also studied. In particular, a new feature regarding 
negaton-positon (negaton-soliton) interaction is analyzed. 

The paper will be organized as follows. We first review the Toda lattice hierarchy in 
section 2. In section 3, through the nth-constrained flow for Toda lattice hierarchy we give 
the Toda lattice hierarchy with self-consistent sources and its Lax representation. Based on 
the Darboux transformation for Toda lattice equation, we develop a method to construct 
the FDT and GFDT with arbitrary functions of time for TLSCS in section 4. In section 5, 
we construct some solutions for TLSCS by using GFDT with arbitrary functions of time. 
Some properties of these solutions are analyzed therein. 

2 Toda Lattice Hierarchy 

We first review the Toda lattice hierarchy. Assume / = f(n,t) for n£Z and tel. Define 
shift operator E and difference operator A as follows 



We often denote E k 



(Ef)(n,t) = f(n + l,t), 

(A/)(n,f) = (E-l)f(n,t) = f(n + l,t)-f(n,t). 
k f(n,t) by f (h) . The inverse of E and A are defined as 



(E- 1 f)(n,t) = f{n-l,t), 




Consider the following discrete isospectral problem [17, 19] 



Lip = {vip) (1) + pip + ifj ( 1} = 



(1) 
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where v — v(n,t), p — p(n, t), ip = ip(n,X,t). This equation can be rewritten as the following 
2x2 matrix eigenvalue problem 

=U(v,p,\)V, (2) 

where 

u ^- x) ==(-"■' x- P ) ■ *=(f 

First consider the following stationary zero- curvature equation for generating function T [20], 

T^U-UT = 0, (3) 

where 

+oo +oo / , \ 

i=0 i=0 V Cl " l / 

and ai,bi,Ci are functions of n and t. The first few of these read 

a = ^, 6 = 0, c = 0, 

ai = 0, bi — —1, ci = v ( \ 

(i) , (i) (i) 

a 2 = v, b 2 = -p , c 2 =pv , 

/ . (IK (1) , / (I) 2 . (1) , (2)n (1)/ 2 , , W\ 



In general 

Aa i+1 = p Aai -Ci-bt v , 
o i+ i = p Oj - Aa,, c, = -v b i . 

Define the modification matrix 

A n := diag(6 n+1 + 5,5), 
where 5 is an arbitrary constant. Let 

n 

v n := (A re r) + + A n = r,A n - J + A n , 

i=0 

and 

-% n = V n V, (n > 1). (4) 



Then the compatibility condition of (2) and (4) gives rise to the following zero-curvature 
representation of the Toda lattice hierarchy 

U tn = UV n -V™U, (n>l). (5) 
Hamiltonian form for Toda lattice hierarchy is given by 

with the density H n := — & n+2 /(n + 1) and Hamilton operator 



J :-- 



v (E^ 1 — 1) 
(l-E)v 



For n — 1 in (5) we have the well-known Toda lattice equation 



(-i) 

v t — vp — vp, 

(1) 

Pt = v -v . 



(6) 



Set v := exp (x ( ^ — x^j, p := x t , equations (6) can be represented as 

x tt = exp (x ( 1} — x*j — exp (^x — x W j . (7) 

3 The Toda lattice hierarchy with self-consistent sources 

First we briefly review the nth- constrained flow for Toda lattice hierarchy [15] which is 
defined as the following system 



5H n \ _n_ ( *Aj 
Sv 
8H„ 



where 



+ E s% =0 (8a) 

Sp / j=i \ Sp / 

L(p j+ := {v(j) j+ ) W + pcj) j+ + (j)^ = \j4>j+, (8b) 

L*4>j- ■= v<j> { p + p4>j_ + <j>f_ = j — 1, . . . ,N, (8c) 



(SXj/Sv, 5\j/5p) T = <t>j-<t>j+\ 
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Define 

n 

f »=E 



j=0 



Ci &i 



) A " i + g^(A^)(lI<' 



/ (1) , (1) 
<t>j-<t> j+ 

'j- ( Pj+ v ( t ) j- ( i ) j+ 



The Lax representation for (8a) is [15] 



(9) 



Following the idea in the continuous case that by treating the constrained flows of soliton 
equation as the stationary equations of the soliton equation with self-consistent sources, we 
define Toda lattice hierarchy with iV self-consistent sources as the following system 



1 

p t 


= J 


Sv 
SH„ 
Sp 


L(j>j+ 


= Xj, 






= x jt 





sr^N SXj 
+ 2^j=l Sp 



J 



(-1) / . a'-- 1 ' 1 A 



U n+1 ' Z^7 = l 



N 

j=l <Pj-Vj+ 



(10a) 

(10b) 
(10c) 



j = l,...,N. 

The Lax representation for (10a) can be obtained from the Lax representation (9) 

= UV, -% n = (A n f n + A n ) tt, (11) 

where A n = diag(6 n+ i — J^jLi + ^' with 5 an arbitrary constant. 

For n = 1, we get the Toda lattice equation with TV self-consistent sources (TLSCS) 



v t = v\p 



/ N \ / N 

= 4i+E^-W)-« (1) (i+E^. 

V .7 = 1 ' ^ .7 = 1 



(1) 



(12a) 
(12b) 



= Aj-0 j+ , (12c) 
LVi- = A^_, j = l,...,AT. (I2d) 

The Lax representation for TLSCS is obtained from (11) by taking n = 1, and <5 = A/2. 
We prefer to rewrite it equivalently in the following scalar form 

Xip (13a) 



= v ip + p - ^ + W 



N 



Lcj) j+ = \j<j> j+ , 



N. 



(13b) 

(13c) 
(13d) 
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With the substitutions 

N 

v := exp(x X) - x), p ■= x t - ^ 4>j+4>j-> ( 14 ) 

3=1 

equations (12a) and (12b) can be written as 

N 

( (-1) \ / - . ' 

X tt 



= exp (x ( 1} - + fa+ti- ) 

V j= i J 

(IV \ i\ 

1 + Yl ^jlfa- ) + Yl 

.7=1 7 7=1 



AT \ TV 

(^■-) t . (15) 

3=1 7 J=l 



4 The Darboux Transformations for TLSCS 

4.1 The forward DT with arbitrary functions of time 

Based on the Darboux transformation for Toda equation in [17, 19], we can find the following 
theorem. 

Theorem 1. Given the solution v, p, x, <pi± (i = 1,...,N) for (12) and (15), and eigen- 
function for (13a) and (13b), let f and g be two independent eigenf unctions of (13a) 
and (13b) with X — /i. Denote h = f + a(t)g with the coefficient a(t) being an arbitrary 
differ entiable function oft. Then the FDT is defined as follows 

m=^-^\ (16a) 



h m 

o 

~h 2 



, (i) , (-i) 

v[l]=v wh ^—, (16b) 



h h ( - n 



P[±]=P-j-m+—> (16c) 



1 , ( // 



2 °V/i ( 

h (J: 



x [l]= x " +-\og( — ) , (16d) 



i+ [l] = i+ - — <j> i+ , (16e) 



A- 1 E(/i0 t _) + ^(t) . 

&-[l] = ^ , i = l,...,N (16f) 

W, + [l]=c(t)(/-^/ (1) ), (16g) 

W+i,-[l] =d(t)\, (16h) 
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where c(t) and d(t) are arbitrary fixed differ entiable functions of t satisfying c(t) ■ d(t) = 
-a/a, Ki (t) = ^[v {1) h m cf>i- - ^_]U=o- Namely v[l\,p[l\, x[l],&±[l] (i = 1, . . . , N + 1), 
Av+i = fJ> and ip[i\ give a new solution to (12) or (15) and (13) with N + 1 self- consistent 
sources. 

This theorem can be proved by straightforward calculation. 

Remark 1. Theorem 1 serves as a non-auto-Bdcklund transformation between two TLSCSs 
with degrees of sources N and N + 1 

We give an example for obtaining solution via Theorem 1. 

Example 1 (Rational solutions). Starting from trivial solution v = l,p = and vanishing 
sources, the Lax pair (13) becomes 



Lip = ifj W + ip^ X) 



(i) 



(17a) 
(17b) 



Let ip = exp(kn — e k t) be the solution of (17) w.r.t. A = 2cosh(/c) ; / = (an — at + b) exp(— t), 
g = exp(— t) be independent solutions of (17) w.r.t. A = 2, where a6l- {0} and b G K are 
arbitrary constants. Let h — f + a(t)g with the differentiate function a(t). Then we obtain 
the rational solution with a pair of non-vanishing sources 



1-e 



MA 
v[l] 

P[l] 
x[l] 

1+ [1] = -c(t) 



k an + b — at + a(t) 



exp(kn — e k t), 

an + a + b — at + a(t) J 
[an + a + b — at + a(t)} [an — a + b — at + ot(t)\ 

[an + b — at + a(t)] 2 
an — a + b — at + a{t) an + b — at + a(t) 



an + b — at + a(t) an + a + b — at + a(t) ' 



^log 



an + b — at + a{t) 
an + a + b — at + a(t) 
ae~ l a{t) 
an + a + b — at + a(t) ' 



0i- [i] = d(ty 



an + b — at + a{t) ' 
where c(t) and d(t) satisfying c(t)d(t) = —a/a. 
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4.2 The Multi-time Repeated FDT 

Theorem 2 (The multi-time repeated FDT). Given the solution v, p, x, <pi± (i = 
1,...,N) for (12) and (15), and eigenfunction ip for (13a) and (13b), let fj and gj be 
independent eigenf unctions of (13a) and (13b) w.r.t. distinct fij (j = Let ctj(t) 

(j = 1, . . . , I) be arbitrary smooth functions of t. Denote hj = fj + ctjgj- Then the l-times 
repeated FDT is given as 



p[i 



X 

<f>i- 

<pN+j - 



cas (ip,h u ...,h i) 
cas (hi, ...,hi) 



(i) 



to cas (hi, . . . , hif ) cas (hi, . . . , h{) ' 



cas (h^..., hi)' 



p + 



cas (hi, ...,hi) 



(-i) 



cas (hi,..., hi) 



cas (hi,..., hi) cas(/ii,...,/i,) (1) 
cas (hi, ... , hi) 



(0.1, 
x f - lo_ 

1 [cas (hi, . . . , hi) 
cas((j) i+ ,h 1 ,. . . ,hi) 

cas (hi, . . . , hif 5 
cas(0j_, hi,..., hi) 
cas (hi, . . . , hi) 
cas (fj, h,..., hi) 



(i) 



i = l,...,N 



Cj(t) 



cas (hi, . . . ,h t ) 



(i) ' 



dj(t)- 



cas (hi, . . . , hj, . . . , hi 
cas (hi, ...,h t ) 



(i) 



3 = 1,. ..,1 



where 



cas (hi,..., hi) :=det yh) i= i i > 

/ j=o,..'.,i'-i 

is £/ie Casorati determinant, and 

cas (hi, ...,h t ):= det . =1 ^ 



j=0,2,...,Z 



cas (<f>i-, hi, ... , hi) := det 



Si(hi) 
hi 



(i) 



(i) 



*1 



(i-i) 



(i-i) 



(18a) 
(18b) 
(18c) 

(18d) 
(18e) 
(18f) 
(18g) 

(18h) 
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where 



Si(hj) := S((f)i-, hj) = A 1 E((p i ^h j ) + (fij - Xi) 1 [v W hf (p^ - hj(f)fl}\ n=0 . 

The hj means the removal of this term and Cj(t) ■ dj(t) = (—iydj/atj. Functions v[l], 
p[l], x[l], <f>i±[l] (i = 1, . . . , N + I) and \N+j = (j = 1, • • • , satisfy (13), (12) and (15) 
with N replaced by N + I. 

Sketch of Proof . For the proof of formulae (18a), (18b), (18c), (18e), (18g) see [19]. The 
formula (18d) is proved by using (14). Formulae (18f) and (18h) are proved by induction. 
The calculation is lengthy but rather straight forward. We omit it. 

□ 

Remark 2. The multi-times repeated FDT (18) provides a non-auto-Bdcklund transforma- 
tion between two TLSCSs with degrees N and N + I. 

4.3 The generalized forward Darboux transformations 

It is well known from [17] that the positon solution of the Toda lattice equations are obtained 
by computing the limit ki — > k\ in the result of two-step DT, where k\^ are parameters 
of eigenfunction with which DT is generated. In order to construct positon solutions for 
TLSCS, similar consideration can be made in our case of FDT with arbitrary functions of 
time. However, in our case, the arbitrary time functions aj(t) must be carefully chosen to 
balance the divergence of sources. Possible candidates for otj{t) are exponential functions. 
By using it, we get the following GFDT. 

Theorem 3 (GFDT). Given the solution v, p, x, 4>i± (i — 1, . . . , N) for (12) and (15), and 
eigenfunction ip for (13a) and (13b), let Fj, Gj (J — 1, ... ,1) be pairs of independent eigen- 
functions of (13a) and (13b) corresponding to distinct X^+j. Let f r , g r (r = 1, . . . , I) be pairs 
of independent eigenfunctions of (13a) and (13b) corresponding to X^+i+r = ^N+i+ri^r) , 
where X^+i+r(^r) is analytic function of parameter u r e C. Let (Xj(t) (j = 1, . . . ,/), f3 r {t) 
(r = 1, . . . , J) be arbitrary differ entiable functions of t. Let m r G N, m r > 2 (r = 1, . . . , I), 
m := (mi, . . . , mi). Denote qj = Fj + a>j{t)Gj, h r = f r + g r respectively. Then the following 
GFDT with I times of FDT with arbitrary time functions aj(t) and I times of generalized 
FDT of multiplicities m r with arbitrary time functions f3 r (t) is given by 

#,m] = %M (iga) 
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X 



<t>i- 

<PN+j,+ 



I, m 
p[l, m 

/, m 
I, m 

I, m 
I, m 



(/>N+j- 
4>N+l+r,+ 

(pN+l+r- 



l, m 
I, m 

/, m 



= v 



A (1) A (_1) 



A? 



'(-i) 



A/, m Aj 



(1) • 

m 



(I+|m|) 1 / A; m 

x + 2 log Ur 



Aj, m (0. 



i+) 



A, (1) ' 
A«, m (0i_) 



A, 



i = l,...,iV 



= 9(0 



A, m (^) 



Af 



(1) ' 

H,ra 

(i) 



A^ m 

A,, m (/ r ) 
" AS, ' 

Pr A , 
A,m 



J = l,...,/ 



r = l,...,J 



where Cj(t) ■ dj(t) = (—lydj/ctj, |m| := X] m r- Symbols are defined as follows: 
Az, m (» :=cas(^,gi,g 2 , 



A, m : = 



cas(gi,g 2 , 

.,c / " 1 ^+(- i r / " 1 ^/), 



hi,d u}I h I ,dl i h I , . 



Ai >m : = cas(gi,g 2 , 



•,^'- 1 fc/ + (-l) m ' -1 )9/(*)^), 



(19b) 
(19c) 

(19d) 

(19e) 

(19f) 
(19g) 

(19h) 
(19i) 

(19j) 



Af m := cas(gi,...,qj,...,gi; 
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h, a^h, a^h u . . . , o^-'h + {-ir^mgu • • • ; 

hi, d^hj, d 2 hj, . . . , d^- l h: + {-l) mi -%(t) 9l ) (1 < j < I), 



A \% '■= cas(gi,g 2 , 

h, d Ul h u & ui h u . . . , d^h + (-l)^- 1 p 1 (t)g 1 ; • • • ; 
h 8 h B 2 h B mr ~ 2 h 

hj, d Ul h T , dlhj, . . . , d^hi + {-l) mi - x fr(t) gi ) (1 < r < I), 
where cas and cas are defined in Theorem 2. 



.) := det ^Wj(gi), . . . , Wi(gi); 



Wi(/n), ^Wi^i), ^w^/ii), . . . , ^T 1 1 - 1 w i (/n) + (-^^-^iWi^i); • • • ; 

w^/ij), c^w^), c£ /Wi (/i 7 ), . . . , d^r^iihi) + (-l) m7 " 1 /3/w i (^)) , 

where for any solution ip of (13a) iw^/i eigenvalue \, Wi(tp) is I + |m| dimensional column 
vector defined as 

Wi(^) = (^(^),^ (1) ,...,V <i+|m '" 1) ) T ) 
and «s a scalar defined as 



SM) := S(0,_, V) = A- 1 ^^) + (A - A^V V'V*- - 



n=0- 



TTien m], v[l, m], m], x[Z, m], 0j±[/, m] and Aj (i — 1, . . . , N + l + I) satisfy (13), (12) 
and (15) wift AT replaced by N + I + I . 

Proof. Without loss of generality, we only prove the special case 1 — 0,1 — 1. The multiplic- 
ity mi is denoted by m (m > 2). And \\(uji) is denoted by A (a;), which is analytic function 
of parameter oo. f\, gi, hi and 0i(t) are denoted by /, g, h and /3(t) respectively. 

Let uj s = uj + ee s , where e s (s — 1, . . . , m) are distinct complex constants, e is a small 
parameter. Define g s (t) = exp(Q s b s (t)) , where b s (t) are arbitrary differentiable functions of 
t satisfying Y17=i b s(t) = f3(t), and 

tt s = - — -t If (Ui-Us), p s = - ( TT7 IT ( e i~^)- 

(m-1)! (m-1)! 
We have the obvious relation fi s = e m ~ 1 p s . 
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We have the following important observations. Let u(A), v(A) be two column vectors 
of dimension m, u(A), v(A) be two column vectors of dimension m — 1, whose components 
are analytic functions of A = X(uo). Denote u s = u(\(lu s )), v s = v(A(u; s )), u s = u(A(cj s )), 
v s = v(A(cj s )) then 



det(ui + gi vi, . . . , u m + QmV m ) 

-i) , 

II ( e i ~ e *) det ( u + v,d w (u + v), 



^.|m(m— 1) 



1! • • • (m 

l<i<j<m 



. . . , d™-\u + v) + (-l) m - l f3(t)^ + o(^ m(m - 1} ), (21) 

det(Ui + 01 Vi, . . . , Ur^To^Vr, . . . , u m + fav m ) 
£ i(m-2)(m-l) 



H (ej - e t ) det (u + v, d u (u + v), . . . , «9^ 2 (u + v)) 



l!...(m-2)! 1 

+ ( £ |(m-2)(m-l)^ (22) 

These observations can be easily obtained by insert the Taylor expansion of u s , v s and 
v s at e = into the determinant. 

By applying multi-time repeated FDT (Theorem 2) with h(X(u s )) = f(X(u s ))+Q s (t)g(X(u s )), 



c s = —1, d s = (— l) s_1 0s/0s ( s — l)---; m )> we obtain V[m], t>[m], p[m], x[m], (pi±[m] 
(i = l,...,N + m) and Atv +s = X(cu s ) satisfying (13), (12) and (15) with N + m self- 
consistent sources. Since such solution contain an arbitrary parameter e, letting e — > will 
also give rise to solution to (13), (12) and (15) with N + m self-consistent sources. So by 
doing this with use of (21) and (22), we get the special case (I = 0, m = m) of formulae 
(19a), (19b), (19c), (19d), (19e) and (19f). We also obtain m self-consistent sources 

cas(/,/i,...,3T 1 >0 



<pN+sA m \ = 



cas(/ i ,...,,9™- 1 /i+(-l)— ip(t)g)' 



cas(/i,...,c>-- 2 /i) 

Qn+s -I'm] = S — — ; r T7TT\ n i s = 771. 

caa(h,...,d^- 1 h + (-l) m - 1 /3(t)gY ' ' 

Since 07v +Si+ [m]'s are equal, 07v +Si _[m] differ only in coefficients b s and they all corresponding 
to eigenvalue X(cu), it is reasonable to combine m self-consistent sources to one self-consistent 
source by denoting 

0JV+i,+ [O,m] = <f> N+li+ [m], 

m 

<f> N+1 -[Q,m] = ^2(f> N+s _[m]. 

8 = 1 
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Thus we arrive at (19i) and (19j). 

The general case can be proved by repeating above procedure. 



□ 



5 Solutions of TLSCS 

The GFDT technique enable us to construct various types of solutions to the TLSCS. Starting 
from trivial solution v — 1, p — for (12) with iV = 0, by choosing specific solutions of Lax 
pair (17) and specific /, I and m in Theorem 3, we can construct multi-solitons solutions, 
(multi-)positon solutions (or higher order), (multi-)negaton solutions (or higher order) and 
(multi-)soliton-positon solutions, (multi-)soliton-negaton solutions, (multi-)positon-negaton 
solutions etc.. 



5.1 Soliton solutions 



Let F = z/™exp(n7 — ise 1 t),G = z/ ra exp(— wy — z/e~ 7 t) be solutions of (17) with respect 
to A = 2z/cosh(7), where 7 G R\{0}. Let a = exp(— 2a(t)), where a(t) is an arbitrary 
differentiable functions of t, v — ±1 . Define 

q = F + aG = 2u n exp(— z/cosh(7)t — a) cosh([/), U := wy — v sinh(7)t + a. 

Then using theorem 3 with I — 1, I — 0, we have the 1-soliton solution for (12) or (15) with 
N = 1. 

Ak cosh(f/) 



^ so1 = exp(ikn - e t) 



1 - ve l 



cosh(C/ + 7) 



v so1 = 



sol 



cosh(C/ — 7) cosh([7 + 7) 
cosh 2 ([7) ' 
cosh(C/ — 7) cosh(C/) 



sol 



X 



0+' 



01°' = 



V- 



log 



— z/- 



cosh(C/) cosh(C/ + 7)' 

cosh(C/) 
cosh(C/ + 7) 

exp(— z/ cosh(7)t — a) 



-z/ n sinh(7)c(t) 
v n d(t) 



cosh(/7 + 7) 
exp(z/ cosh (7) t + a) 



2 cosh(/7) 
where c(t)d(t) = 2d. 

The 1-soliton solution of TLSCS is similar to the 1-soliton of ordinary Toda lattice 
equation [21]. The main difference between them is that the travelling speed of 1-soliton 
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solution for TLSCS can vary with t, for the travelling speed is i/sm1 ^ 7 a . The multi-soliton 
solutions for TLSCS (12) can be constructed similarly. We omit it. 

5.2 Positon solutions and negaton solutions 
5.2.1 One-positon solution 

Let/ = ci exp(— t cosu>) cos(nu — t sinu;), g = c<i exp(— t cos u) sin(nu; — t sinu;) be solutions 
of (17) w.r.t A = 2cosu;. Set constants C\ = cos 9, c 2 = sin9, where 9 ^ ^tt {n G Z). Define 

h = f + g = exp(— t cosu>) cos(Y), Y = nu — t sinu; — 9, 

and 

Z = d u Y = n-tcosiu, rj = Z + ~p(t) sin(20). 

2 

where f3(t) is an arbitrary differentiable function. Then using theorem 3 with the specification 
I — 0, / = 1, mi — 2, one gets one-positon solution for TLSCS (12) with N — 1. 

yPOS = [(V + |) sin^ + | sin(2r + 3a;)] [(y - \) sing; + \ sin(2F - uQ] 

[(y + \) sinu; + \ sin(2F + uo)] 2 
P os = 77sin(2a;) + sin(2y) _ (y + 1) sin(2uQ + sin(2r + 2uQ 
P ' (y + |) sinu; + |sin(2F + u;) (y + §) sinu; + \ sin(2F + 3a;) 1 ' 



x pos = - log 



I r (y + i) sin u + \ sin(2F + u) 



-i 2 



(77 + f ) sinu; + \ sin(2F + 3a;) 



(23c) 



V+ V ; (??+§) sinu; + ±sin(2y + 3u;) V ; 

= _^ exp(tcosu;)cos(y + c) 

(77 + |) sinu; + |sin(2F + u;) v ; 

The scattering properties can be analyzed, resembling with [17]. It is noticed that the 
scattering matrix is an identity matrix, which is the basic feature for positon solutions. 

Then we discuss the analytic properties of one-positon solution. An investigation on 
positon profile x pos shows that the one-positon profile decays to zero slowly at ±00, oscillates 
and possesses singularities during their propagation. The singularities of one positon profile 
are determined by zeros of 

A pos := (77 + -) sinu; + - sin(2F + uj). 

The speed of positon profile are defined by the speed of singularities [16]. Assuming 
n (t) to be one of the places where singularity occurs, differentiating by t on the two side 
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of the equationA pos (n , t) = 0, one finds that the singularity propagation is governed by a 
nonlinear ODE 



sinu 



n = 



UJ 



1 + - 



uj cos uj — \uj(5 sin(2#) — sin a; 



sin uj + uj cos(2n uJ — 2t sin uj — 29 + uj) 



The singularities of + os and 0?° s are determined by zeros of (AP os y and Ap° s , respec- 
tively. 

5.2.2 Two-positon solution 

Let fi := cos 9 i exp(—t cos uJi) cos(nuji — tsinuji), ^ := sin6'jexp(— tcosu;,) sin(nu;; — tsinc^), 
% — 1,2 be solutions of (17) w.r.t A; = 2coscjj, 9i ^ (n G Z). Define 

hi = fi + 9i = exp(-t cos^) cos(Yj), Y, L = nuj,i - tsincj; - 9 h 

Zi = d w Yi = n - t cos uji, ^ = Z { + ^Pi(t) sin(26>;). 

Then using theorem 3 with / = 0, I = 2, m = (2, 2), one obtains the 2-positon solution for 
TLSCS (12) , 



x 



2p 



log 



cas (cos Yi, r/i sin Yi, cos Y 2 , r\ 2 sin Y 2 ) 

cas (cos Yi, 771 sin Yi, cos Y 2 , r\ 2 sin Y 2 ) 

cas (sin Yi, cos Yi, rji sin Yi, cos Y 2 , 772 sin Y 2 ) 



0i+ = --sin(20i)exp(-tcoswi) 

1 cas (cos Yi , 771 sin Yi , cos Y 2 , r^ 2 sin Y 2 ) 



0^. = —Pi exp(t cos c<jj 



cas (cos Yi, cos Y 2 , 772 sin Y 2 



(i) 



cas (cos Yi, rji sin Yi, cos Y 2 , 77 2 sin Y 2 ) ' 

and 02± have the similar formulae. x 2p describes a wave profile oscillates, decreases to zero 
slowly as |n| — > 00. X + and behave like 0{n~ l ) as |n| — > 00. 

The Positon-positon interaction can be analyzed as follows. Fixing 771, assume j^l — ► 00 
when t — > ±00, then 



x 2p ~ - log 



1 



^fa + |) + ^sin(2Yi + 3m) 
A(r)i + |) + B sin(2Y! + 5wi) 



, as £ — > ±00, 



where 



A := 2 sinu;i sinu;2 + sin(3u;i) sinc^2 — 2 sin(2a;i) sin(2o; 2 ) + sincui sin(3a; 2 ), 

3 1 
B := sinuj 2 cos(2c<Ji) + - sincj 2 — sin(2w2) cosc^i + - sin(3c^2)- 
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Thus, we observe one-positon in terms of rji and Y\ at t — ±00. And there is no phase shift 
during the interaction. Analogous, fixing r\ 2 we observe one-positon at t — ±00 in terms of r/ 2 
and Y 2 . Similarly, there is no phase shift in the course of interaction. Positons propagations 
transparently as if others were absent. This transparency of interaction is a remarkable 
feature of positon solutions. The N-positon solution for TLSCS (12) is obtained by taking 
fi := cos 6i exp(— t coscui) cos(ncUi — t sin a;,), g t := sin0jexp(— tcosc^) sin(nu;j — tsincuj), 
% — 1, . . . , N be solutions of (17) w.r.t A, = 2 cosu;j, 0j 7^ ^ir (n G Z), = + gi and use 
theorem 3 with I = 0, I = N, m = (2, 2, • • • , 2) e N*. 

We can construct N-negaton solutions to TLSCS (12) by considering eigenfunctions /; : = 
e"exp(nt<;j — e^e^t + 0,), ^ := e™exp(— nwj — e;e~ Wl t — 9 { ) of (17) with A; = 2e i coshc<; i , 
distinct ^ 6 1, ^, € R, = ±l,hi = fi + gi and using theorem 3 with I = 0, I = N, 
m= (2,2, ••• ,2) G N N . 

The GFDT is quite general to construct not only N-soliton (positon, negaton) solutions, 
but also solutions of combined type. 

5.3 Posit on- negaton solution and interaction 

Let 

fi := cos^i exp(— t coscjx) cos(no;i — tsinc^x), 
g 1 := sin^x exp(— t cosux) sin(no;i — tsmui), 
f 2 := e n 2 exp(nu 2 - e 2 e w H + 9 2 ), 
g 2 := e 2 exp(-nw 2 - e 2 e~ W2 t - 9 2 ), 

where 9i, 9 2 , uji and uj 2 are real numbers, 9\ 7^ kir/2 for any integer k. Define 

hi — fi + g 1 — exp(— t cosc^x) cosYi, h 2 = f 2 + g 2 = 2e^ exp(— te 2 cosh uj 2 ) cosh Y 2 , 
where 

Yi := nui — t sin Ui — 9 1 , Y 2 := nu 2 — te 2 sinh u 2 + 9 2 . 

Denote 

Vi ■= d^Yi + -(3i(t) sin 20x, ?? 2 := d U2 Y 2 + -(3 2 (t). 

Then according to theorem 3 with the specification / = 0, / = 2, m = (2, 2), we obtain the 
positon-negaton solution for TLSCS (12). 
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To analyze the interaction, we always assume fait) and fait) tending to =1=00 as n — > ±00. 
And all parameters, including oj\ uj 2 and e 2 are positive. Fixing 771, if for particular fait), r\ 2 
increases faster than sinh(2Y" 2 ) (i-e. \rj 2 / sinh(2F 2 )| — > 00) as t — > 00, which indicates that 
negaton travels at a very high speed dominated by fa(t) then 

\ sin(2Y"i + 3cji) + sin u;i (771 - A 8 ) 



x 



pm 



log 



1 - As 



where 



_ I sin(2Yi + 5cji) + sin u x (771 

cos a*! — 2e 2 cosh w 2 — 2 cosh w 2 cos uj 1 + 2 



ast — >• ±00, 



e tJ2 (e 2 e- W2 - e iWl )(e 2 e~ UJ2 - e^ 1 ) 
Thus we see the one-positon profile (see (23c)) with neither phase shift nor displacement at 
the two end. That is to say the negaton is transparent for positon, which is a phenomenon 
never observed in the ordinary Toda lattice case. 

If r/ 2 increases slower than sinh(2F 2 ) (i.e. |r/ 2 /sinh(2Y" 2 )| — > 00) as \t\ increases, which 
implies a slowly traveling negaton in comparison with the previous case, then 

2 

- 2u; 2 , 



x 



pn 



1 



~ - log 

2 5 



sina;i(?7i + A 9 ) + \ sin(2Y"i A u 



x 



pn 



~ - log 



sin wi (771 + 1 + A 9 ) + !sin(2Yi + 3cji - A n ) 

sin ^1(771 + Ai ) + \ sin(2Yl + 5wi + A n ) 
sincj!^! + 1 + A10) + |sin(2Yl + 1uj x + A n ) 



as t 



-00, 



-, 2 



+ 2u 2 , as t — > +00, 



where 



-e 
2 



-2uj 2 



no 



D 2u> 2 



(e 2 e^ 2 - Ze^)ie 2 e^ - Ze~^) - 4 
(e 2 e^ 2 - e iaJl )(e 2 e~ aJ2 - e^ 1 ) ' 
(e 2 e-^ 2 - 3e iwi )(e 2 e- W2 - 3e~ iu!1 ) - 4 
(e 2 e W2 - e iu)1 )(e 2 e^ - e~ iu}1 ) '' 



and 



^11 •- 



« log 



e 2 e^ 2 - e iuJl 



e 2 e" 2 - e~ lu)1 

are all real constants. Thus positon travels with phase shift determined by A 9 , A 10 and An, 
with displacement determined by 4cj 2 in the course of collision. This is general phenomenon 
caused by existence of negaton. 

If we fix a coordinate frame which travelling with negaton profile, then 

-1 2 



x 



pn 



log 



I sinh(2Y" 2 + 3u 2 ) + (r/ 2 + 3) sinhw 2 
\ sinh(2Y" 2 + 5uj 2 ) + (i] 2 + 4) sinhu; 2 



as t — > ±00. 



Thus, the negaton travels insensitive about the existence of positon. 
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Soliton-positon and soliton-negaton solutions can obtained in the same way by using 
theorem 3 with I — 1, I — 1, m = 2, q :— F + aG, h := fi+gi and q := F + aG, h := f2 + Q2 
respectively. <7j i — 1,2 are defined in the beginning of this section and F, G, a are defined 
in the section 5.1.) 

6 Conclusions 

Based on the constrained flows of Toda lattice hierarchy, we constructed Toda lattice hier- 
archy with self-consistent sources and their Lax representation. 

We developed a method to construct FDT with arbitrary functions of time and the 
GFDT with arbitrary functions of time, which, in contrast with the well-known Darboux 
transformation for Toda lattice, provide non-auto-Backlund transformation between two 
TLSCSs with different degrees and enable us to obtain various explicit solutions to TLSCS. 
Resembling the ordinary Toda lattice case, this system possesses solutions of rich families, 
including solitons, positons, negatons and the solutions of combined types. A number of 
solutions are listed by our method. The investigation on these solutions shows a quite 
similar nature with solutions of ordinary Toda lattice. However, the new feature concerning 
interactions between negaton and positon (or soliton) which is different from the ordinary 
Toda lattice case is also stated. This difference is caused by the wide range of variation of 
the speed of negatons in TLSCS cases. We note that variation of speed is a common feature 
for continuous and discrete systems with self-consistent sources, see [6, 7] etc. 

It is convinced that our approach for constructing systems with self-consistent sources 
and generalized forward Darboux transformation technique are available for other discrete 
systems. Some investigation will be present in the forthcoming paper. 
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